Matter Creation via Vacuum Fluctuations in the Early Universe and Observed 

Ultra-High Energy Cosmic Ray Events 
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Cosmic rays of the highest energy, above the Greisen-Zatsepin-Kuzmin cut-off of the spectrum, may 
originate in decays of superheavy long-living X-particles. These particles may be produced in the 
early Universe from vacuum fluctuations during or after inflation and may constitute a considerable 
fraction of Cold Dark Matter. We calculate numerically their abundance for a wide range of models. 
X-particles are considered to be either bosons or fermions. Particles that are several times heavier 
than inflaton, mjnflaton w 10 13 GeV, and were produced by this mechanism, can account for the 
critical mass in the Universe naturally. In some cases induced isocurvature density fluctuations can 
leave an imprint in anisotropy of cosmic microwave background radiation. 

PACS numbers: 98.70.Sa, 95.35.+d, 98.80.Cq 



I. INTRODUCTION 

According to the Greisen- Zatsepin-Kuzmin § (GZK) 
observation, the energy spectrum of Ultra High En- 
ergy (UHE) cosmic rays produced at extragalactic dis- 
tances should exhibit an exponential cut-off at energy 
E ~ 5 x 10 10 GeV. However, a number of cosmic ray 
events with energies well beyond the predicted GZK cut- 
off were recently observed by the various experimental 
groups |2j. This is an obvious contradiction with the 
standard cosmological and particle physics models and 
clearly requires some new physics beyond the Standard 
Model. 

Conceptually, the simplest explanation |^,|| could be 
that the highest energy cosmic rays are produced in the 
decays of heavy long- living particles in cosmologically lo- 
cal part of the Universe. We will call these progenitors of 
the cosmic rays as X-particles. The mass of AT-particles 
has to be very large, mx ^ 10 13 GeV, for them to be re- 
sponsible for the cosmic rays events in the energy range 
E > 10 11 GeV. It was noticed [|]|] that sucn heavy 
particles are produced in the early Universe from the 
vacuum fluctuations and their abundance can be correct 
naturally, if the standard Friedmann epoch in the Uni- 
verse evolution was preceded by the inflationary stage. 
Basically, it is the same process which during inflation 
had generated primordial large scale density perturba- 
tions and seeded formation of galaxies and galaxy clus- 
ters. However, there is a difference. Inflationary stage 
is not required to produce superheavy particles from the 
vacuum, unlike to the case of generation of the long wave- 
length perturbations. Rather, the inflation provides a cut 
off in excessive production of heavy particles which would 
happen in the Friedmann Universe if it would start from 
the initial singularity ||. Note that we consider sterile 
or almost sterile X-particles, so that their production in 



usual plasma interactions and decays can be neglected. 
Thermal production of heavy X-particles was discussed 
in Refs. pj M a nd production during reheating considered 
recently inR]. 

Particle production from vacuum fluctuations in the 
Friedmann Universe during matter or radiation domi- 
nated stages had been considered long ago |^,[| and this 
gives the basis for simple estimates of X-particle abun- 
dance JBj which have the mass of oder or smaller than 
the Hubble constant at the end of inflation, mx < Hi. 
However, the most interesting case of heavier X-particles, 
mx > Hi, requires detailed calculations. Relevant calcu- 
lations were done already by Chung, Kolb and Riotto || . 
However, to describe the end of inflation and the transi- 
tion to the Friedmann universe, they match the fixed de 
Sitter background to the subsequent radiation or matter 
dominated expansion either as an instantaneous junction 
or with the help of a smoothing function. The result is es- 
pecially sensitive to the details of the junction procedure 
in the case of large mx which is of interest for us here. 
Further, the range of models considered in Ref. |l was 
restricted to the case of scalar particles with conformal 
coupling to gravity. 

The purpose of the present paper is to calculate pro- 
duction of superheavy particles from the vacuum in the 
inflationary Universe for a wide range of models, and to 
give detailed and extended discussion of our other re- 
sults already reported in letter j5|]. Considering parti- 
cle creation, we do not make any approximations. To 
this end we find numerically the exact evolution of the 
scale factor in the model of the "chaotic" inflation p| . 
Our basic formalism relies on the method of Bogolyubov 
transformations, which for the case of particle creation by 
non-stationary gravitational field was developed in refs. 
pO|,p|. We consider the possibility for X-particles to be 
either bosons or fermions. 

In Section II we remind the reader about the explicit 
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relations between mode functions of the field, coefficients 
of the Bogolyubov transformation and the particle num- 
ber. In Section III we present results of our numerical 
calculations of particle production from the vacuum fluc- 
tuations both for Friedmann and inflationary stages in 
the universe evolution. In Section IV we relate results of 
previous section to the present-day density of X-particles. 
In Section V we discuss relevance of our results to the 
UHE cosmic rays events, and Section VI contains our 
conclusions. 



II. PARTICLE PRODUCTION BY A 
NON-STATIONARY GRAVITATIONAL FIELD 

Here we summarize the basic formalism of gravita- 
tional particle creation in an expanding universe which 
we employed in our analysis. For more details see Refs. 
|js|,fio| . We choose the metric to be conformally flat at an 
early cosmological epoch, ds 2 = a(n) 2 (di] 2 — dx 2 ). The 
number density of Wparticles created from the vacuum 
in a time varying cosmological background can be written 
as 
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where j3k are the Bogolyubov's coefficients which relate 
"in" and "out" mode functions, k is the comoving mo- 
mentum, and is the sum over spin states. The ex- 
pression ([!]) gives the number density of particles only, 
with the equal amount of antiparticles being created in 
the case of charged fields. The creation of Bose and Fermi 
particles is to be considered separately. 

i) Bosons. The mode functions, \k = Xfcl 7 ?) °f a scalar 
Bose field are solutions of the oscillator equation 



x" k + ^kXk = o , 
with the time dependent frequency 

uj 2 = k 2 -— (1-60+ ™W 
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Here ' = d/drj. The constant £ describes direct coupling 
to the curvature, with £ = corresponding to the mini- 
mal coupling and £ = 1/6 being the case of conformal 
coupling. Equations for massless conformally coupled 
quanta are reduced to the equations in flat space-time 
and such particles are not created. For massive particles 
conformal invariance is broken and particles are created 
regardless of value of £. Given the initial (vacuum) con- 
ditions 



(4) 



the Bogolyubov's coefficients at any time moment n are 
found to be 



2 |x' fc | 2 +c 2 | Xfc | 2 -2c 
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ii) Fermions. The relevant mode functions of the Fermi 
field satisfy the oscillator equation with the complex fre- 
quency 



x"k + ~ im x a')xk = 



(6) 



where the real part of the frequency is given by uj 2 = 



k 2 + m\a? . We choose 



Xfc (0) = Wl-^, X , k (0) = -i"Xk, (7) 
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as the initial conditions. In this case we find per spin 
state 



i fl i2 - m x a - Im(xfcXfc 
^ 2^ 



(8) 



Summation over spins gives the factor of 2 in Eq. ([!]). 



III. NUMERICAL RESULTS 

A. Power-Law Cosmology 

In this section we consider particle creation during the 
regular, i.e. non-inflationary, stage of the Universe ex- 
pansion. We shall refer to this stage as the Friedmann 
stage of the Universe evolution. During the Friedmann 
stage the scale factor is given by the expression a(t) cx t a 
with a < 1. This case is of interest since a = 1/2 
and a = 2/3 correspond to the radiation and matter 
dominated stages of the Universe expansion, which are 
inevitable. Moreover, in some cases (and we shall en- 
counter such situations below) particle creation during 
this stage gives a dominant contribution, while creation 
during inflation (or at the transition between stages) is 
negligible. 

Let us first consider massive particles conformally cou- 
pled to gravity. It is the particle mass which couples 
the system to the background expansion and serves as 
the source of particle creation. Therefore we expect 
nx oc m x a~ 3 at late times when particle creation di- 
minishes and the number of particles conserves in a co- 
moving volume. Because the scale factor has an arbi- 
trary normalization, for definiteness it is more convenient 
to rewrite this expression using relations a oc (mt) a , or 
a oc (m/H) a . We prefer the latter choice since the exact 
value of time depends on the previous history in the infla- 
tionary cosmology and consequently the time coordinate, 
t, is ambiguously defined. Moreover, it is a non-zero value 
of the Hubble constant which is responsible for particle 
creation in an expanding universe. The particle produc- 
tion is expected to stop when H mx- Therefore, at 
the late epoch we can parameterize the anticipated for- 
mulae for nx as 
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where the constant C a depends only upon the back- 
ground cosmology, i.e. a, and it can be found numer- 
ically. Note that for the radiation dominated case it 
was found in Ref. f§, n x « 5.3 x 10- 4 m 3 x (m x t)-^ 2 . 
We confirm this result and extend it here to a broad 
range of cosmologically interesting values of a, which in- 
cludes the matter dominated case as well. Our results 
are summarized in Fig |l| where we plot C a as a func- 
tion of a. The particle number reaches this asymptotic 
value, Eq. (||), if the initial vacuum state was defined 
at the sufficiently early epoch, H(0) ^> m X - We found 
that -ff(O) 10 4 mx is sufficient for the results to be- 
come practically independent of the chose of the initial 
time if a — 1/2, with smaller (larger) ratio of H(0)/mx 
required at larger (smaller) a. 
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FIG. 1. Coefficient C' a , denned in Eq. (9), is shown as 
a function of a for background cosmology with a power law 
scale factor a oc t a . 

For the radiation dominated Universe one finds fix = 
Px/pc — ™>xnx 327rGi 2 /3 with the present value of fix 
being equal to fix ~2x I0^ 2 (m x /Mp l )^mxt e , where 
t e is the time of equal densities of radiation and matter in 
the ft = 1 Universe. This gives fl x ~ (m x /10 9 GeV) 5 / 2 . 
We see that stable particles with mx 10 9 GeV will 
overdose the Universe even if initially they were in a 
vacuum state and were created from the vacuum during 
the regular Friedmann radiation dominated stage of the 
evolution. (It is possible to separate the vacuum creation 
from the creation in collisions in plasma since X-particles 
may be effectively sterile.) 

However, this restriction will be not valioQ if the evo- 
lution of the Universe, as it is believed, was more com- 
plicated than the simple radiation dominated expansion 
from a singularity. The Hubble constant may have never 
exceeded mx, which is the case of inflation, H(0) w 
rricj, w 10 13 GeV. Moreover, compared to the case con- 
sidered above, the density of X-particles created during 
inflation is additionally diluted by the entropy release in 



'Another possibility is the late time entropy release. 



reheating after inflation. 

B. Inflationary Cosmology 

Any viable modern cosmological model invokes the hy- 
pothesis of inflation (for a review and list of references, 
see refs. [[l2],[l3|). During inflation the Universe expands 
exponentially which solves the horizon and flatness prob- 
lems of standard Big-Bang cosmology. Inflation is gen- 
erally assumed to be driven by the special scalar field 
(j> known as the inflaton. Fluctuations generated at in- 
flationary stage may have the strength and the power 
spectrum suitable for generation of the large-cale struc- 
ture. This fixes the range of parameters of the inflaton 
potential. For example, the mass of the inflaton field has 
to be ~ 10 13 GeV. During inflation, the inflaton field 
slowly rolls down towards the minimum of its potential. 
Inflation ends when the potential energy associated with 
the inflaton field becomes smaller than the kinetic en- 
ergy. Coherent oscillations of the inflaton field contained 
all the energy of the Universe at that time. It is possible 
]l4j that a significant fraction of this energy was released 
to other boson species after only a dozen oscillations of 
the inflaton field, in the regime of a broad parametric 
resonance. This process was studied in detail [ ^5|Jlr| . It 
was shown that even rather heavy particles with masses 
by an order of magnitude larger than the inflaton mass 
can be produced quite copiously (see also |L7| ) . Applying 
these results to the case of our interest, we find that the 
stable very heavy particles, m^ < m x ^5 10 m^, will be 
generally produced in excess in this process. 

However, if the parametric resonance is ineffective for 
some reason, and one estimates the particle number den- 
sity after inflation at the level of initial conditions used in 
Refs. |l^|, one finds that fix might prove to be of the ac- 
ceptable magnitude. This level is saturated by the funda- 
mental process of particle creation during inflation from 
vacuum fluctuations and it is the same process which 
generated primordial large scale density perturbations. 
Parametric resonance for X particles is turned off if X 
field is either a fermion field or its coupling to inflaton is 
small, g 2 < 10 4 (m x / m^) 4 (m^ / M P {) 2 |]|. (Parametric 
creation of fermions was considered in Refs. |l8). While 
it is not a "resonance" process, it might be useful for 
creation of very heavy fermions.) 

Particle creation from vacuum fluctuations during in- 
flation (or in the de Sitter space) was extensively studied 
pj| , pO| , usually in the case of small m x and in applica- 
tion to generation of density fluctuations necessary for 
the large scale structure formation. The characteristic 
quantity which is usually cited, the variance of the field, 
(X 2 ), is defined by an expression similar to Eq. (jl]). 
In the typical case the difference is given by the factor 
2 sin 2 (wfc?7)/ajfe in the integrand. For example, for the 
scalar Bose field with the minimal coupling to the cur- 
vature, (X 2 ) = 3H?/8ir 2 m 2 x if m x < H % For 
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massless self-interacting field (X 2 ) « Q.132Hf/y/\ [|T). 
Particle creation in the settings close to our problem, but 
for the specific case of the Hubble dependent effective 
mass, mx(t) oc H(t), was considered in Ref. p2| . 

In general, if rax ~ -Hj ~ we can again expect 
that the particle density will be given by Eq. (||). In gen- 
eral, the coefficient C in this equation will be a function 
of the ratio Hi/mx, of the coupling constant £, and will 
depend on details of the transition between inflationary 
and matter (or radiation) dominated phases, etc. For 
mx > Hi an exponential suppression is expected. De- 
tailed numerical study of X-particle abundance created 
from the vacuum in inflationary cosmology is the subject 
of this Section of the paper. 

The case of large mx with application to the dark mat- 
ter problem was already considered in Ref. ||. However, 
the analysis was restricted to the case of conformally 
coupled scalar particles only, and approximations were 
made to the extent that the fixed de Sitter background 
was matched to the subsequent radiation or matter dom- 
inated expansion either as an instantaneous transition or 
with the help of some smoothing function. Particle cre- 
ation in the case of large mx can be especially sensitive 
to (some extent arbitrary) details of such junction pro- 
cedure. 

In the present paper, considering particle creation as 
described by Eqs. (|lj)-(|j), we do not make any approxi- 
mations. We assume a specific model of inflation, namely 
"chaotic" inflation pj| , and find numerically the exact 
evolution of the scale factor in this model. We consider 
both the "massive" inflaton with the scalar potential 
V{4>) — m 2 (f) 2 /2, and the "massless" inflaton with the 
potential V(4>) = A0 4 /4. The normalization to the ob- 
servable large scale structure requires m 2 /Mpj 10~ 12 
in the former model and A ~ 10 -13 in the latter model. 
We consider creation of conformally and minimally cou- 
pled scalar X-particles, as well as creation of X-fermions. 

Our results are summarized in Fig. ^ (massive infla- 
ton) and Fig. || (massless inflaton). We define the energy 
density in X-particles at late times as px = rnxnx- In 
the case of the massive inflaton we normalized px by 
the total energy density. The ratio of those two quan- 
tities reaches some asymptotic value and then remains 
constant in a matter dominated universe, and can be 
measured at any sufficiently late time in the numerical 
simulation which does not include the inflaton decay. For 
each value of mx we integrated equations (|j)-@ until 
the asymptotic constant value of this ratio was reached. 
Similarly, the convenient quantity to be measured in the 
massless inflaton model is proportional to px /p 3 ^ 4 , which 
also becomes constant at late times. The curves labeled 
by £ = and £ = 1/6 correspond to the minimal and the 
conformal coupling to gravity, respectively. The dashed 
line describes creation of fermions. 

We see that particle production is most effective in 
the case of scalar particles with minimal coupling, ex- 
ceeding px generated in other models by many orders of 



energy density in X-particles to the total energy den- 
sity is independent upon mx at mx < m in this model. 
This might seem as easily understandable. Indeed, mag- 
nitude of fluctuations generated during inflation in de 



Sitter space is (X 2 ) 



x 



if m x < Hi. Multiply- 



ing this by m x one would find px ~ m x (X ) — const, 
seemingly in agreement with observed behavior, Fig. |^. 
However, this agreement is rather coincidental since it is 
possible to relate p and (X 2 ) in this simple manner at 
late times only. 
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FIG. 2. Ratio of the energy density in X-particles to the 
total energy density at late times in a model with the massive 
inflaton, V{4>) = m 2 </> 2 /2, as a function of X particle mass, 
mx- For the inflaton mass we defined mis = m/10 13 GeV. 
The dotted line is the low mass asymptotic, Eq. (Mj. 
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Note also that the ratio of the 



FIG. 3. Energy density of created X-particles in the 
massless inflaton model, V(4>) = A</> 4 /4. We defined 
A13 = A/10" 13 . 

In our numerical calculations we had found that the 
variance (X 2 ) of the field X measured at the end of in- 
flation is independent upon mx if the mass of X is small. 
At some later epoch when H w mx, which will be long 
after the end of inflation if X is a light field, the field X 
starts to oscillate on all scales, including k = 0. Only at 
this time, which we denote by tx, all field fluctuations are 
transformed into non-zero particle density and we can use 
Px =m x n x ~ m 2 x 



(X 2 ). The variance of X fluctuations 
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was unchanged on large scales, starting from the end of 
inflation down to the time tx- So, when the field starts to 
oscillate, px oc m\. However, the energy density of the 
inflaton field, p = 3H 2 /8ttG decreased during this time 
interval in proportion to H 2 (t x )/H 2 (0) w m x /H 2 (0). 
That is why the ratio of the energy density in A-particles 
to the total energy density does not depends on mx when 
measured at t > tx- 

Variance of the field X is different from the usually 
calculated for the fixed de Sitter inflationary background 
because we consider the actual evolution of the scale fac- 
tor and the value of the Hubble parameter is not con- 
stant during inflation, being larger at earlier times. Cor- 
respondingly, the number of created particles per decade 
of k grows logarithmically towards small k if mx is small. 
(Power spectrum behaves similarly.) For definiteness, 
we restricted our calculations to ~ 10 15 for the expan- 
sion factor during inflation. The examples of the particle 
number, k 3 nx(k) for several values of mx are shown in 
Fig. |^ at the moment corresponding to 10 completed in- 
flaton oscillations. The particle momentum is measured 
in units of the inflaton mass. In contrast to this, in fixed 
de Sitter background Att 2 {X 2 ) « H 2 J d\nk {k/Hf-' 2v 
with 3 — 2v w 2mj c /3H 2 at small mx and consequently 
(X 2 ) cx \/m 2 x . When our code was tested on fixed de 
Sitter background it reproduced the proper power spec- 
trum exactly. Note that the power spectrum in fixed de 
Sitter background grows towards large values of k, which 
is opposite to the behavior of Fig. 0. 
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FIG. 4. Spectrum of created particles, k 3 n(k), in a model 
with massive inflaton is shown for several choices of the mass 
of scalar X-particle with the minimal coupling (solid lines) and 
the conformal coupling (dotted line). Masses and momenta, 
k, are given in units of the inflaton mass. 

Therefore, calculations which would be based on cus- 
tomary procedure of matching a fixed de Sitter back- 
ground to a subsequent Friedmann stage would give 
wrong results in this case. Note, however, that the re- 
sult depends upon duration of inflationary stage, so our 
results in the case of scalars with minimal coupling and 
mx < m should be considered only as a lower bound on 
X-particle abundance. 



Matching is also dangerous in the case of large mx, 
and its influence was studied in detail in Ref. ||. When 
the change is too abrupt, it generates artificial particles. 
This can easily happen for mx > m, see e.g. |23|] where 
excessive production was found. At mx ^ m the number 
of created particles decreases exponentially with mx ■ For 
the case of massive inflaton our result shown in Fig. || can 
be fitted as 

?j « 10" 10 ^TO? 3 exp(-27rmx/m) , (10) 

where m^ = m/10 13 GeV. (Note that the fitting formu- 
lae is approximate and arbitrary to some extent). 

As Fig. H shows, the power spectrum of fluctuations 
in X-particles is almost scale independent at small k if 
mx/m ~ 1. Therefore, if such particles constitute a con- 
siderable fraction of dark matter, this fluctuations will 
be transformed into isocurvature density perturbations 
at late times and can affect large scale structure forma- 
tion. Isocurvature fluctuations produce 6 times larger 
angular temperature fluctuations in cosmic microwave 
background radiation (CMBR) for the same amplitude of 
long-wavelength density perturbations compared to the 
adiabatic case p4| |. Therefore, to fit observations by a 
single spectrum, the mass fluctuation spectrum in isocur- 
vature cold dark matter cosmology must be tilted (with 
respect to scale invariant spectrum) to favor smaller 
scales. 

Let the power spectrum of the field fluctuations be 
k 3 Px (fc) oc k@ . To fit to the second moments of the large- 
scale mass and (not contradict to) cosmic microwave 
background distributions requires (3 > 0.25. Models with 
(3 ranging from 0.3 to 0.6 were constructed and considered 
in Refs. [g5|,|6). It is interesting that in our case the spec- 
trum is indeed correctly tilted, see Fig. ^, for mx which 
is somewhat larger than the inflaton mass, On the small- 
est scale of our integration, k w 10 -15 we find (3 w 0.1 
for mx/m = 2 and (3 ~ 0.3 for mx/m — 3. Since (3 
decreases with the length scale increasing, on scales cor- 
responding to the current horizon (3 will be smaller. Be- 
cause of that the amplitude of perturbations will be also 
too small on cosmological scales when (3 will be reaching 
desired magnitude, (3 > 0.25. Indeed, density fluctua- 
tions, Spx/px, which are roughly proportional to spec- 
tra shown in Fig. 0^ are reaching unity near the break 
in the power spectrum (fc ~ 1 at the end of inflation) 
and on the scales of interest amplitude of fluctuations 
became too small already at mx/m = 3. However, their 
magnitude may be just right, Spx/px ~ 10~ 5 , on scales 
k > 10~ 25 for the case mx/m ss 2 (or slightly larger) 
if fix ~ 1, see Fig. ||. Magnitude can not correspond 
to the observable on the galaxy scales since then there 
will be too much fluctuations in CMBR. This gives a 
constraint on the models. X-particles with minimal cou- 
pling to gravity, mx > 2m and Qx ~ 1 are excluded. 
(Such particles with small contribution to O may exist 
though.) However, the magnitude of the density fluctu- 
ations induced in the process of X particle creation can 
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correspond to the observable on the horizon scale, and be 
responsible for fluctuations in CMBR. X-particles with 
mx /m ~ 2 and Qx ~ 1 can give such contribution natu- 
rally. (Source for the galaxies seeds has to be different in 
this case.) Note that the density autocorrelation function 
is proportional to the square of the field autocorrelation 
function in the present case, and the density fluctuations 
are non-Gaussian. 

In the case of fermions or scalar particles with confor- 
mal coupling, the number density of created particles has 
a kink at mx = mm the model with massive inflaton, see 
Fig. ||. This signifies that different mechanisms are re- 
sponsible for particle creation at mx > m and mx < m. 
Particles heavier than the inflaton are created during in- 
flation, while particles lighter than the inflaton are cre- 
ated after inflation during the regular Friedmann stage of 
expansion. We see that particle creation during this lat- 
ter stage is actually more effective (nx rapidly increases 
when mx drops slightly below m). Hence, we can assume 
that the number density of created particles at mx <C m 
is given by Eq. (| with C a = C 2/3 . (C 2/3 w 9.7 x 1CT 4 
for the conformal scalars, see Fig. El) Since the total 

" — 3H /8nG, we find for 



energy density is given by p 
a = 2/3 



p x _ 8?r m\ 



3/2 



(11) 



This function is plotted in Fig. ||by the dotted line which 
shows that already at mx < 0.5 m the number density 
of created particles is given by Eq. (||), indeed. Note 
that in Ref. || all particle creation was attributed to 
the transition from the inflationary to matter dominated 
phase, while we find that for mx < m the dominant 
particle creation occurs at a later stage, when H w k w 
mx (for mx < m). The correct interpretation allows us 
to identify the low mass asymptotic, Eq. (O), see also 



Scalar particles with minimal coupling to gravity are 
produced predominantly during the inflationary stage for 
all mx, and the resulting number density is much higher, 
see Fig. || 

Conformal coupling, £ = 1/6, is a stable point of renor- 
malization group equations. Otherwise £ is a running 
coupling constant and there is no reason to expect the 
coupling to be, say, minimal at mass scales of inflation. 
The dependence of X-particle abundance on the coupling 
to the curvature scalar is shown in Fig. |5[ 

The situation in the case of the massless inflaton, 
V{4>) — A</> 4 /4, is quite similar. However, now there is 
no distinctive mass scale in the equations of motion for 
the scale factor, and the kink in the function nx(mx) 
does not appear. The quantity \/A0 might be consid- 
ered as an effective inflaton mass, but it changes during 
the evolution. For definiteness we plot nx as a func- 
tion of mx / "v^Mpi- Note that the transition from the 
inflationary to the radiation dominated stage occurs at 
(f> ~ 0.35Mpi 1TH] and this may define the more natural 



mass scale for the ratio mx/m^n ((/)). Indeed, we observe 
that nx reaches its maximum when mx ~ 0.35\/AMpi, 
see Fig. |[ 

We find that the number density of produced fermions 
in both cases is approximately equal to the number den- 
sity of scalar particles with conformal coupling. The main 
difference amounts to a spin factor of two. 
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FIG. 5. Ratio of the energy density in X-particles to the 
total energy density at late times in a model with the massive 
inflaton, V((f>) = m 2 <f) 2 /2, as a function of £. 



IV. PRESENT DENSITY OF X-PARTICLES 

Let us estimate the present day number density of X- 
particles. Consideration is somewhat different for mas- 
sive and massless inflaton models. 
1. Massive inflaton, V(<fi) — m|0 2 /2. 
In this case the inflation is followed by the stage of mat- 
ter domination which lasts until inflaton decays and ther- 
mal equilibrium is established. Let us denote the temper- 
ature at the beginning of radiation dominated stage as 
T*. Actually, a transition to the radiation dominated uni- 
verse does not necessarily occur with the final state be- 
ing in thermal equilibrium. Complications arise if there 
are light bosons in a theory, raj <C m^, even relatively 
weakly coupled to the inflaton, g 2 > 10 4 m|/M|i - kH. 
In this case the inflaton will decay via parametric reso- 
nance after just a several oscillations and relativistic de- 
grees of freedom start to dominate the equation of state. 
This happens typically when the energy density in the 
inflaton oscillations is red-shifte d by a factor r w ICP 6 



compared to the value m^Mpj jl4 15 1. Matter is still 
far from being in thermal equilibrium, but since Uni- 
verse expansion already at this point changes to radia- 
tion dominated it is convenient to characterize this radi- 
ation dominated stage by an equivalent "temperature", 
T* ~ r 1 / 4 ^/m^A/pi. In any case, T* is the model depen- 
dent parameter and we leave it as such. 

Until the moment when the equation of state changes 
from matter dominated to radiation dominated, the ratio 
of the energy density in X-particles to the total energy 
density retains its value reached at the end of inflation, 
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px/p = const. This ratio is plotted in Fig. g for different 
models as a function of model parameters and it defines 
the present day ratio of px to the critical density in the 
Universe, p c , as follows: 



fi x h 2 = 22L 



P -fT 



•fink 2 



(12) 



where T is present temperature of cosmic microwave 
background, and flnh 2 = p ra dh 2 /p c « 4.31 x 10~ 5 is the 
fraction of critical density that is in radiation today [|~3| . 
The following definitions are used here p c = SHqM^/Stt 
and Hq = lOO/i km sec -1 Mpc -1 . Parameter 7 describes 
the change in entropy per comoving volume for the whole 
cosmological history starting from T = T* down to the 
present day, i.e. so = 7 s* • In simple cosmological mod- 
els, 7 ~ 1. However, in models with subsequent latent 
heat release (due to phase transitions or due to decay of 
long living particles, e.g. moduli fields) 7 can be very 
large. 

Let us assume for illustration purposes that T» = 10 9 
GeV, i.e. (T*/T )(Q R h 2 ) » 1.8 x 10 17 . (This is the 
highest allowed reheating temperature when gravitinos 
in supergravity theories are not overproduced.) With 
7 ~ 1 wc find that fixh 2 ~ 1 if mx/rn w 2 in models 
with scalar X-particles with conformal coupling, or if X- 
particles are fermions, see Fig. |^. For minimally coupled 
scalars this value of fix is reached for mx/rn ps 3. 

To summarize this subsection: fixh 2 is found 
multiplying the data of Fig. || by a factor 1.8 x 
10 17 (T„/ 7 10 9 GeV)m 2 3 . 

1. Massless inflaton, V(<f>) = A<f> /4. 

Before and after the decay of the inflaton the Universe 
expands as radiation dominated. In Section III for the 
state which precedes the decay of inflaton oscillations we 
had calculated the quantity px/VAMpip 3 ^, where p is 
the inflaton energy density. After the decay of the infla- 
ton, in a state of thermal equilibrium we can use 
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1/4 



P 



3/4 



(13) 



where is the number of relativistic degrees of freedom 
at temperature T*. Assuming instantaneous transition 
between those states, = p, and again using fix = 
{px/T s )(T s /p c ) together with s = 7 s * > we nnd 



fl x h 2 



pi 



30 



-1/4 



V\Mp\p 3/4 



n R h 2 . 



(14) 



We see that VI xh 2 can be found by multiplying the data 
plotted in Fig. || by a factor 2.5 x 10 20 A^/7<72ooi where 
we defined 3200 = g*/200, and used V~\Mpi/To w 1.6 x 
lO 25 ^/^. Note that y/\M vx « 3.9 x 10 12 VA~ 13 GeV, 
so that already with 7 = 1 we shall have fix < 1 if 
m x > 4 x 10 13 GeV. 



V. X-PARTICLES AND UHE COSMIC RAYS. 

Why were the Ultra High Energy cosmic rays observed 
above the Greisen-Zatsepin-Kuzmin energy cut-off? One 
solution to the problem might be provided, for example, 
by the existence of some exotic particles which are able 
to propagate (evading the GZK bound) from cosmologi- 
cal distances and yet interact in the Earth's atmosphere 
like a hadron. A particle with such conflicting properties 
was found in a class of supersymmetric theories |^7| . Al- 
ternatively, high energy cosmic rays may have been pro- 
duced locally within the GZK distance. One possibility 
is connected with the destruction of topological defects 
p8[ , while another one is connected with decays of pri- 
mordial heavy long- living particles The candidate 
V-particle must obviously obey constraints on mass, den- 
sity and lifetime. 

In order to produce cosmic rays in the energy range 
E > 10 11 GeV, the mass of V-particles has to be very 
large, m x > 10 13 GeV [|||]. The lifetime, Tx, cannot be 
much smaller than the age of the Universe, r f=a 10 10 yr. 
With such a short lifetime, the observed flux of UHE 
cosmic rays will be generated with the rather low density 
of V-particles, fix ~ 10 -12 , where fix = mxnx/Pcxiti 
nx is the number density of X-particles and p cl -i t is the 
critical density. On the other hand, X-particles must not 
overdose the Universe, fix ^ L With fix ~ 1, the 
X-particles may play the role of cold dark matter and 
the observed flux of UHE cosmic rays can be matched 
if tx ~ 10 22 yr. The allowed windows are quite wide 
||, but on the exotic side, which may give rise to some 
problems. 

The possibility that X-particles were generated from 
the vacuum by gravitational interactions |^|Jg| may an- 
swer to the question why their abundance should lie in 
the interesting range fix ^ 1 and was considered in detail 
in the present paper. 

The problem of the particle physics mechanism respon- 
sible for a long but finite lifetime of very heavy particles 
can be solved in several ways. For example, otherwise 
conserved quantum number carried by X-particles ma 
be broken very weakly due to instanton transitions 
or quantum gravity (wormhole) effects [Q. If instantons 
are responsible for A"-particle decays, the lifetime is es- 
timated as tx ~ rn^} ■ exp(47r/ax), where ax is the 
coupling constant of the relevant gauge interaction. The 
lifetime will fit the allowed window if the coupling con- 
stant (at the scale mx) is ax « 0.1 0. 

A class of natural candidates for superheavy long- living 
particles which arise in string and M theory was re- 
evaluated recently in Ref. |3(J and particles with desired 
mass and long life-time were identified. Another inter- 
esting candidates were found among adjoint messengers 
in gauge mediated supergravity models |j 31 1 . 
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VI. CONCLUSIONS 

We have shown that the very weakly interacting su- 
perheavy X-particles with mx = (a few) • 10 13 GeV may 
naturally constitute a considerable fraction of Cold Dark 
Matter. These particles are produced in the early Uni- 
verse from vacuum fluctuations during or after inflation. 
Related density fluctuations may have left an imprint in 
fluctuations of cosmic microwave background radiation 
if scalar X-particles with minimal coupling to gravity 
are approximately twice heavier than the inflaton and 
tlx ~ f ■ Decays of X-particles may explain UHE cosmic 
rays phenomenon. 

Our hypothesis has unique observational consequences. 
If UHE cosmic rays are indeed due to the decay of these 
superheavy particles, there has to be a new sharp cut-off 
in the cosmic ray spectrum at energies somewhat smaller 
than mx- Since the number density nx depends ex- 
ponentially upon mx/nij,, the position of this cut-off is 
fixed and can be predicted to be near w 10 13 GeV, 
the very shape of the cosmic ray spectrum beyond the 
GZK cut-off being of quite generic form following from 
the QCD quark/gluon fragmentation. Next generation 
experiments, like the Pierre Auger Project [|32f| . High 
Resolution Fly's Eye Q, the Japanese Telescope Ar- 
ray Project [34|| , may prove to be able to discover this 
fundamental phenomenon. 

We conclude that the observations of Ultra High En- 
ergy cosmic rays can probe the spectrum of elementary 
particles in the superheavy range and can give an unique 
opportunity for investigation of the earliest epoch of evo- 
lution of the Universe, starting with the amplification of 
vacuum fluctuations during inflation through fine details 
of gravitational interaction and down to the physics of 
reheating. 
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